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Abstract 



A test on the numerical accuracy of the semiclassical approximation as a 
function of the principal quantum number has been performed for the Pullen- 
Edmonds model, a two-dimensional, non-integrable, scaling invariant per- 
turbation of the resonant harmonic oscillator. A perturbative interpretation 
is obtained of the recently observed phenomenon of the accuracy decrease on 
the approximation of individual energy levels at the increase of the principal 
quantum number. Moreover, the accuracy provided by the semiclassical ap- 
proximation formula is on the average the same as that provided by quantum 
perturbation theory. 
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Recently, there has been considerable renewed interest in the various as- 
pects of the Semi-Classical Approximation (SCA), a powerful motivation be- 
hind that being the problem of the so-called quantum chaos (see for example 
references [1,2,3,4,5]). An important aspect is represented by the quantum 
energy levels, and in this connection one recent work [6] shows that the pre- 
dictions of individual levels by SCA (by this we mean the Bohr-Sommerfeld 
formula, or one of its generalizations to the non-integrable case, such as 
EBK; see e.g. [3,4]), worsen as the quantum number increases, contrary to 
the naive expectation. We argue that this result can be interpreted as follows: 
if h, no matter how small, is kept fixed, the SCA on the individual levels has 
the meaning of a perturbation theory (PT) in h. Therefore the accuracy of 
the approximation decreases for higher levels (to get good agreement it is 
necessary, as is well known, to implement the classical limit h — > 0, n — > oo, 
nh — I classical action; see e.g. [7]). 

The aim of this paper is to clarify this point, from the theoretical stand- 
point and from the computational one as well, considering a scaling invariant 
potential, which makes ordinary quantum PT strictly equivalent to a power 
expansion in h. We do actually observe that, for H fixed, the perturbation 
strength has to be decreased to keep the accuracy at a constant value as 
the quantum number increases; however we also observe that the algorithm 
provided by the appropriate SCA is always comparable to the algorithm pro- 
vided by ordinary quantum PT. A very good agreement between the "exact" 
eigenvalues, obtained by numerical diagonalization of the Schroedinger op- 
erator, and the semi-classical ones, is indeed observed in presence of high 
unperturbed degeneracy. 
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The most significant examples to carry out this comparison are repre- 
sented by non-separable two dimensional systems exhibiting both regular 
and irregular spectrum [8], i.e. in particular, non-uniform behaviour of the 
level spacing, and among these the simplest one is the Pullen-Edmonds model 
[9]. Its quantum hamiltonian is: 

where m has been put equal to 1. For x = (1) reduces to a resonant two- 
dimensional harmonic oscillator of levels {n\ + n 2 + l)h = mih, mi = 1, 2, 
of multiplicity mi. 

The scaling transformation qi — > Vhqi, q 2 — > Vhq 2 yields the unitary 
equivalent operator: 

H = - ^) + \{ql + ql) + xHql]. (2) 

The coupling constant has become xh, which clearly shows equivalence be- 
tween expansions in powers of x or of /i (an analogous result holds for any 
other polynomial perturbation). Moreover, the symmetry of the potential 
enables us to split the hamiltonian matrix, computed on the harmonic os- 
cillator basis, into submatrices reducing the computer storage required. The 
matrix elements of (1) can be written: 

< n x n 2 \E\mn 2 >= h{ni + n 2 + 1)5 „' ini <W 2 + - ^) 5 n[ ni -2 + 

+ ^ni + l){n l + 2)5 <ni+2 + (2m + l)S n[ J x (3) 
*[\fn 2 (n 2 - l)5 n ' 2 n 2 - 2 + \j(n 2 + l)(n 2 + 2)<^„ 2+2 + (2n 2 + 1)^„ 2 ] 

and each submatrix can be labelled by the parity of the occupation numbers 
ni, n 2 . We restrict from now on to the invariant subspace spanned by mi 
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even, i.e. rt\ and n 2 of opposite parity. The eigenvalues of H in this subspace 
have constant multiplicity 2 [9]. Therefore the level m\h = 2sh splits into s 
levels for x > 0. 

The appropriate SCA is here provided by the Bohr-Sommerfeld quanti- 
zation of the resonant (or secular) canonical perturbation theory [10], also 
known, in this particular the Birkoff-Gustafson normal form [11,12], 

which we now construct at first order. Starting from the classical Pullen- 
Edmonds hamiltonian: 

Hd = \(pI+ v\) + \{q\ + ql) + xqhl (4) 

we introduce the standard action-angle variables (/, 0) by the canonical 
transformation: 

qi = ^JWi cos 0i 



i = 1,2. (5) 
Then (4) becomes: 

H d = h + I 2 + A X lili cos 2 6 1 cos 2 2 . (6) 

The second canonical transformation into the well known "slow" and "fast" 
variables: 

'a, = i 1 + i, U = * + * (7) 

A 2 = h-h { 02 = 01 - 02, 

eliminates the dependence on the "slow action" A 2 in the unperturbed part, 
so that the hamiltonian becomes: 

H d = A l + X (Ai - Aj) cos 2 (0! + 2 ) cos 2 (0x - 2 ). (8) 
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We now eliminate the dependence on the angles up to order x 2 by resonant (or 
secular) canonical perturbation theory [10]. To eliminate the dependence on 
the "fast angle" 0i it is enough to average the perturbation on this variable. 
This yields: 



2tt 

and thus: 



1 r 2?r 1 

-/ d0 x cos 2 (0 X + 2 ) cos 2 (0x - 2 ) = -(2 + cos40 2 ), (9) 

!7T JO o 



H d = A X + |(A 2 - A 2 2 ){2 + cos40 2 ). (10) 



The dependence of 2 on the perturbation part can now eliminated by a 
further canonical transformation. The Hamilton- Jacobi equation for the 
perturbation part is in fact: 

[A 2 -(— ) 2 ](2 + cos40 2 ) = ^ (11) 

0(p2 



83 ±, 



A 2 (2 + cos40 2 ) -K 



d4>2 \ 2 + cos 40 2 
and thus the Hamiltonian (9) becomes: 



(12) 



H d = B 1 + ^K(B 1 ,B 2 ), (13) 

where: 

1 r (99 

Bl = Au B2 = _/^ 2 _ (14) 

It appears from the structure of equation (12) that the motions generated by 
the perturbation part of our system have the following qualitative behaviour: 

< K < B\ rotational motion 
K — B\ separatrix (15) 

B\ < K < 3Bf librational motion. 
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The appearance of a separatrix (which is not immediately obvious in the 
(jp,q) coordinates) accounts as is well known (see e.g. [3]) for the stochastic 
layers originating near it. This corresponds to local irregular behaviour of 
the quantum spectrum; one of its manifestations is (see D. Delande in [4]) 
the local shrinking of the level spacing and the tendency to avoided crossing. 
The shrinking of the level spacing is best accounted by the SCA, as we will 
discuss below. 

On the separatrix we have: 



5 1 2 (2 + cos40 2 ) = K, (16) 



while in general: 



2 fb Bj(2 + cosAx)-K 
B * = ± «LH 2 + cos4x - (17) 



where: 



a = 0, b = £ rotational motion 

(18) 

a = (f>-(K, Bi), b = (j>+(K, Bi) librational motion 



<I>±{K,B X ) = ± 7a rccos(-2-2). (19) 



with: 

1 K 
:- arccosl — ^ 
4 B{ 

Now the approximate hamiltonian (13) depends only on the actions so that 
a semiclassical quantization formula for the (mi even part) of spectrum of 
the operator (1) can be obtained by a straightforward application of the 
Bohr-Sommerfeld quantization rules [10]. Set therefore: 

h = (n 1 + l/2)ft (2Q) 
h = (n 2 + l/2)h, 



whence, from (6): 



Set: 



A x = (n 1 +n 2 + l)h 
A 2 = (ni - n 2 )h. 



A\ = mih 
A 2 = m 2 h, 
by comparison of (19) and (20) we obtain: 

mi = ni + n 2 + 1 
m 2 = m — n 2 , 
where mi = 2, 4, ... and m 2 = ±(mi — 1), ±(mi — 3), ±(mi — 5), ... 
Finally: 

i?i = mih, B 2 = m 2 h; 
then the semiclassical approximation to the quantum spectrum is: 

E mi ,m 2 = mih + -K(mih, m 2 h), 

o 

where K is implicitly defined by the relation: 
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I r 

m 2 n = ±— / 

7T J a 



b dx K m ± h ) 2 ( 2 + cosAx ) ~ K 
X \ 2 + cos4x 



and: 



(21) 



(22) 



(23) 



(24) 



(25) 



(26) 



(27) 



a = 0, b = f 0<K < (mih) 2 

a = <j>-{K,B 1 ), b = (f) + (K,Bi) ( mi h) 2 < K < 3(m^) s 
Remark that for \m 2 \ < [ami] we obtain the quantization of the rota- 
tional motions, while for \m 2 \ > [ami] ([x]=integer part of x) we have the 
quantization of the librational ones. Here, by (17): 

2 ff , 11 + cos 4a; 



a 
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dx\ 



2 + cos Ax 



0.602. 



(28) 



Moreover, it immediate to see that for mi fixed the function K, and hence 
the semiclassical energy E mi>m2 , is a decreasing function of the secondary 
quantum number m 2 . It is furthermore proved in [13] that (25) coincides 
with the exact quantum spectrum up to terms of order H and % 2 . The 
numerical computations (see Fig. 4 below) show that at order 1 in x the 
corrections of order % affect at most the eight decimal figure. 

The "exact" levels have been computed, and compared with the semi- 
classical ones as well as with the levels computed by degenerate first order 
quantum perturbation theory [14], for m\ = 1, ...,60 at h = 0.1 and for dif- 
ferent values of x (given the degeneracy, this is equivalent to compute 1800 
different levels). The results obtained for mi = 30, h — 0.1, x — 10~ 3 an d 
mi — 60, h — 0.1, x — 10~ 5 are shown in Figure 1 and Figure 2, respec- 
tively. The local shrinking of the spacing, reproduced by both methods, can 
be immediately noticed; remark that the corresponding semiclassical levels 
are those near the separatrix (by (28), m 2 ~ 18 and m 2 ~ 36, respectively). 

In Figure 3 the function: 

A = \E Ex — E Sc \ (29) 

vs mi is plotted for mix = 1; this shows that, if the coupling constant is de- 
creased in inverse proportion to the principal quantum number the accuracy 
of SCA not only remains constant but actually improves, as anticipated be- 
cause the scaling invariance makes the limit mi ^oo, x ^ 0, m\x^ const 
equivalent to the classical limit m x — > oo, % — > 0, m-Ji — > const. 

In Figure 4 the accuracies obtained thorough semiclassical and quantum 
first order perturbation theories are compared for mi = 60, x — 10~ 5 , ft = 
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0.1, and in Fig. 5 the difference between the two perturbation theories is 
plotted (remark that the energy decreases as m 2 increases). As can be seen 
the agreement with the "exact" levels is very good and the accuracy is on the 
average the same. Remark however that, as it should be expected (the Bohr- 
Sommerfeld rules take no account of tunneling [15]), the lowest accuracy of 
the SCA is reached near vri2 = 36 which corresponds to the levels near the 
separatrix: for those levels the quantum PT is therefore better than SCA. 

The authors are greatly indebted to Dr. Stefano Isola for many useful 
discussions and to Mr. G. Salmaso for his valuable computational assistance. 
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Figure Captions 



Figure 1: Comparison between the "exact" levels (a), the semi-classical ones 
(b), and the levels obtained by first order perturbation theory (c), for \ = 
1(T 3 , h = 0.1, mi = 30. 

Figure 2: Comparison between the "exact" levels (a), the semi-classical ones 
(b), and the levels obtained by first order perturbation theory (c), for m x = 
60, x = 10~ 5 . 

Figure 3: The difference A between the "exact" levels and the semi-classical 
ones vs mi, with m 2 = m 1 — 1 and = 1. 

Figure 4: (a) The difference between the "exact" levels and the semi-classical 
ones; (b) the difference between the "exact" levels and the first order quan- 
tum PT ones; (mi = 60, h = 0.1, \ = f0~ 5 ). 

Figure 5: The difference between the semi-classical levels and the first order 
quantum PT ones; (mi = 60, h — 0.1, x — f0~ 5 ). 
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